Astronomy & Astrophysics manuscript no. cambresy 


©ESO 2011 


January 7, 201 1 





Variation of the extinction law in thie Trifid nebula 

L. Cambresy\ J. Rho^-^, D. J. Marshall^ and W. T. Reach^'^ 

' Observatoire astronomique de Strasbourg, Universite de Strasbourg, CNRS, UMR 7550, 1 1 rue de I'Universite, F-67000 Strasbourg, 
France 

e-mail: cambresy@astro . unistra . fr 
- Infrared Processing and Analysis Center, California Institute of Technology, Pasadena, CA 91 125, USA 

^ Universite de Toulouse, UPS, CESR, 9 avenue du colonel Roche, 31028 Toulouse Cedex 4, France e-mail: 
Douglas . MarshallOcesr . fr 

■* Stratospheric Observatory for Infrared Astronomy, Universities Space Research Association, Mail Stop 211-3, Moffett Field, CA 
94035, USA e-mail: jrho@sofia.usra.edu, wreach@sofia.usra.edu 

Received ; accepted 

ABSTRACT 

Context. In the past few years, the extinction law has been measured in the infrared wavelengths for various molecular clouds and 
different laws have been obtained. 

Aims. In this paper we seek variations of the extinction law within the Trifid nebula region. Such variations would demonstrate a 
local dust evolution linked to variation of the environment parameters such as the density or the interstellar radiation field. 
Methods. The extinction values, A,^IAy, are obtained using the 2MASS, UKIDSS and Spitzer/GLIMPSE surveys. The technique is 
to inter-calibrate color-excess maps from different wavelengths to derive the extinction law and to map the extinction in the Trifid 
region. 

Results. We measured the extinction law at 3.6, 4.5, and 5.8 yum and we found a transition at Ay ~ 20 mag. Below this threshold 
the extinction law is as expected from models for Ry = 5.5 whereas above 20 mag of visual extinction, it is flatter.Using these 
results the color-excess maps are converted into a composite extinction map of the Trifid nebula at a spatial resolution of 1 arcmin. A 
tridimensional analysis along the line-of-sight allowed us to estimate a distance of 2.7 ± 0.5 kpc for the Trifid. The comparison of the 
extinction with the 1.25 mm emission suggests the millimeter emissivity is enhanced in the dense condensations of the cloud. 
Conclusions. Our results suggest a dust transition at large extinction which has not been reported so far and dust emissivity variations. 

Key words. ISM: dust, extinction - ISM: individual object: Trifid nebula - Infrared: ISM 



1. Introduction 

During the last three decades, dust evolution in the interstel- 
lar medium has been studied mostly through the grain emis- 
sion in the far-infrared and then the sub-millimeter wavelengths. 
■ The analysis of the surface brightnesses at 60 and 100 fj.m ob- 
' served by IRAS lead to distinguish a small and a big grain 
component where the small grains are the major contributor at 
60 yum (Boulanger et al. 1988; Laureiis et al. 1991). The transi- 
tion from small to big grains occurs at a typical visual extinction 
of 1 mag. More recently a new transition was discovered in the 
dust composition with the formation of fluffy grains . The analy- 
, sis of PRONAOS sub-millimeter observations by B ernard et alj 
1 d 19991) yielded to the discovery of a population of enhanced- 
emissivity dust grains which have a lower equilibrium temper- 
ature. This emi ssivity enhan cement in the far-infrared by a fac- 
tor of about 3 jCambresv et al. 2001) requires the existence of 
composite grains formed by the coagulation of small grains. A 
large scale study on the whole Galactic anticenter hemisphere 
showed this transition also occurs at low extinction. Ay ~ I mag 
(ICambresv et al.ll2005|). It might be closely related to the recent 
discovery made bv iMolinari et al. (2010a) with Herschel data 
where this same density threshold seems to indicate the transi- 
tion where diffuse clouds start to collapse into filaments. 

With the wealth of available deep data in the near and mid- 
infrared it is now possible to investigate the interstellar grain 
properties at shorter wavelengths. In this spectral domain a lot 



of progresses has been made in the past few years through the 
study of the extinction law. Dense clouds can be probed on 
large scales and for various astrophysical environments. Recent 
studies with Spitzer data (Indebetouw et al. 2005; Flahertv et d] 
'2007';'Roma n-Ztifii^aet al.li2007l: IChapman et a l. 2009) focused 
on different clouds and obtained different extinction laws. We 
propose in this work to look for possible variations of the ex- 
tinction curve within a single cloud. We choose an active star- 
forming region which includes the famous HII region M20, the 
Trifid nebula. Its location in the Galactic plane at only 7 degrees 
from the Galactic center offers a high stellar density which is 
ideal to probe the extinction distribution. The distance of this 
young nebula (10^ years) is still debated. The observation of the 
central 07 star responsible for the photoionization, H D 164492, 
yields a distance from 1.67 kpc dLvnds et al.lll985h to 2.8 kpc 
(iKohoutek et al.|[T999h . This active region has its first genera- 
tion of massive stars interacting with the interstellar medium 
and i s very lik ely triggering a second generation of star forma- 
tion jCernicharo et al. 1998) . The state-of-the-art for the Trifid 
star-fo rming region is presented in a review paper bv lRho et al.l 
(l2008l) . 

The extinction law is generally estimated in molecular 
clouds as a whole and the discrepancies between the different 
results are interpreted as variations of the dust optical prop- 
erties from cloud-to-cloud due to the environment. Obviously 
variations should also be observed within molecular clouds. 
As the density increases from the envelope to the core of a 
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cloud, the dust properties, hence the extinction law, should 
be affected. An attempt to detect such a variation was pro- 
posed by iRom an-Zuniga et al. (2007) in their study of a dark 
cloud core, B59, up to Ay ^ 5 9 mag but they did not suc- 
ceed. lOlofsson & OlofssonI (12010) used five different methods 
to derive the extinction law in B335. The scatter in their re- 
sult was consistent with the uncertainties so they could not con- 
clude on any e v olutio n of the extinction law within the cloud. 
IChapman et al.l (l2009h managed to find variations by averag- 
ing their result from several clouds and they showed evidence 
for grain growth at Ak^ > 1 mag (Ay > 9). A global galac- 
tocentric vari ation of the extinction law has also been reported 
by Zasowski et al.l (l2009h in a large scale study that excluded 
molecular clouds. 

The present study focuses on the variations of the extinc- 
tion law within the molecular cloud associated with the Trifid 
nebul a. Our analy sis relies on the 2M ASS dSkrutskie et alj 
l2006h. UKIDSS (iLawrence et al.1 I200I and GLIMPSE II 
(IChurchwell et al ] |200¥ surveys. Special attention is given to 
the data filtering to limit the possible biases as explained in 
Sect.|2l Our original method to derive the extinction law from 3.6 
to 5.8 yum is based on the color mapping rather than the classical 
analysis of the source list distribution in a color-color diagram. 
The details and the advantages of this technique are presented in 
Sect. [3] while Sect. |4] describes the final extinction map making 
and a tridimensional analysis of the Trifid line-of-sight which 
allows us to disentangle the complexity of this direction. In ad- 
dition we studied the dust emissivity at 1.25 mm by comparing 
the masses derived from the dust continuum emission and the 
dust absorption. We conclude in Sect. |5] 



2. Data preparation 

2.1. Datasets 

Our study covers a large field of 1.6° x 1.4°. The 2MASS neai-- 
infrared data provides JHKs fluxes for the whole sky. The high 
stellar density in the Trifid direction induces a loss of sensitivity 
as 2MASS sources with less than 5" separation are affected by 
confusion. The degradation reaches almost 1.5 mag in each filter, 
with completeness level of only 12.8 mag for the band for 
instance. 

Deeper JHKs fluxes are obtained from the UKIDSS 6th 
data release (UKIDSSDR6plus) of the Galactic Plane sur- 
vey. UKIDSS uses the UKIRT Wide Field Camera (WFCAM; 
Casali et al. '2007). The photometric system is described in 
Hewett et al. (2006), and the calibration is described in 



Table 1. Completeness limit magnitudes. For UKIDSS the lumi- 
nosity function slopes are not well defined, so the peak values of 
the function for sources with an error smaller than 0.15 mag are 
displayed. 



Hodgkin et al.l (I20(j9h . The pipehne processing and science 
archive are described in lHamblyet"an (120081) . 

Longer wavelengths from Spitzer/IRAC are publicly avail- 
able through the GLIMPSE II legacy survey. We used the high 
reliability point source Catalog (v2. 1) which includes data from 
epochs 1 and 2. 

The merged catalog between 2MASS, UKIDSS and 
GLIMPSE tables results from a 1" cone search. Completeness 
limit magnitudes for all filters are presented in Table[T] For bright 
stars, the 2MASS photometry is chosen over the UKIDSS pho- 
tometry which is affected by saturation. Although the UKIDSS 
photometry is directly calibrated on 2MASS we found a system- 
atic difference between the two systems: 

7(2MASS) = 7(UKIDSS) -hO.084 
i/(2MASS) = //(UKIDSS) -0.074 
K,{2M.ASS) = /i:,(UKIDSS) -0.003 
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We corrected the UKIDSS photometry to match the 2MASS sys- 
tem. The UKIDSS astrometry has also been corrected to match 
the GLIMPSE positions by adding 0.12 and 0.07 arcsec to the 
right ascension and declination, respectively. In the following, 
only magnitude measurements with an uncertainty smaller than 
0.15 mag will be considered. It yields about 2 x 10^ sources at 
K,, 6 X 10^ at 3.6 jum in the range R.A. e [269.50,271.29] and 
Dec. € [-23.6,-22.2]. 



2.2. Source selection criteria 

The analysis of the stellar color excess to investigate the cloud 
properties requires all stars to be background stars. The contam- 
ination by foreground objects cannot be neglected for the Trifid 
nebula located between 1.7 and 2.8 kpc probably in the Scutum 
Arm, behind the Sagittarius Arm. The procedure to remove fore- 
ground stars is statistical for two reasons: 1) the number of de- 
tected sources is close to 6 x 10^, and 2) there is a degener- 
acy between intrinsically red foreground objects and blue back- 
ground sources. A color criteria would therefore be inefficient. 
It would actually lead to a source spatial distribution correlated 
with the dust distribution which is not consistent with the sup- 
posedly foreground nature of the stars. The individual identifica- 
tion is generally impossible, except in the high extinction region 
where the degeneracy can be raised. Keeping these constraints 
in mind, we choose a two steps strategy, first purely statistical 
and then individual and manual. 

The statistical step consists in sampling the source catalog 
into a data cube with two spatial dimensions and a third di- 
mension being a color axis. Only one star per box is allowed. 
Basically, stars are sampled by a 4 arcmin grid step according to 
their celestial coordinates and the third axis allows to color-sort 
the stars within a same 16 arcmin^ pixel so that there is only a 
single star per 3D-box. The first slice of the cube represents a 
map where each pixel contains the color of the bluest object of 
the stack. The following slices contain redder and redder color 
values. The completeness limit of UKIDSS being not well de- 
fined we use only 2MASS and GLIMPSE for this statistical pro- 
cedure. The variation of the stellar density makes the number of 
stars per pixel, and therefore the size of the stack, range from 85 
to 387. For a 16 arcmin^ pixel, it means the source density varies 
from 19x10-^ to 87x10"^ deg"^. Once the cube is built, each plane 
from the bluest to the reddest is compared to a raw color map of 
the field, i.e. a non-calibrated extinction map. A correlation with 
the color map starts at the 38* slice, meaning the distribution 
of the ~38 bluest stars in all the 16 arcmin^ pixels is not cor- 
related with the dust distribution. They must be considered as 
foreground objects and the first 38 stars of each pixel, i.e. 8600 
stars deg"^, are thus removed from the original catalog. It repre- 
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sents from 10% to 45% of the sources detected at and [3.6] 
depending on the position, i.e. the local source density. 

At this point, only a fraction of the foreground stars is re- 
moved. A second step is required to filter the catalog in the re- 
gions where the contamination is critical because the foreground 
objects may even be the dominant population. Fortunately, in 
these regions the color degeneracy is raised by the strong red- 
dening from the dust grains. To define these regions we build 
a [3.6]-[4.5] color map following the adaptive method de- 
scribed in Sect. 13.1.11 Then we select the pixels where [3.6]- 
[4.5]>0.35 mag {Ay w 30 mag) which yields a 25 arcmin^ sur- 
face area. Such a high threshold guaranty the degeneracy with 
the intrinsic star color is broken. For this manual operation the 
restriction on UKIDSS data can be waived and ~10^ sources are 
detected in the 25 arcmin area at both and 3.6 with er- 
ror smaller than 0. 15 mag. The spectral energy distribution from 
1.25 jjm to 8 /im is examined for these sources. For each fore- 
ground star candidate, a color-composite JHKs image is built to 
compare the color of the candidate itself with the color of the 
surrounding objects. This visual inspection allows to make the 
final classification of the source since a foreground star appears 
bluer than the surrounding objects. Among the ~10^ sources de- 
tected at Ks and 3.6 /im, 20% are classified as foreground stars. It 
corresponds to a total foreground star density of ~37x 10^ deg"^. 
It worth noting that for sources together detected at J, H and 
instead, the contamination level in the 25 arcmin^ reaches 55%. 
This is because many background stars are lost when the J de- 
tection is required. 

Beside the foreground stars, there is a contamination by the 
Young Stellar Objects (YSOs) that are formed in the molecular 
cloud associated with the Trifid nebula. Several studies, at vari- 
ous wavelengths, provide us with a source list of YSOs that we 
remo ved from o ur catalog. The most important contribution is 
from iRho et all (j2006) who analyzed the infrared colors from 
Spitzer IRAC and MIPS data. It yielded the identification of 161 
YSOs (37 class I/O, 13 hot excess. 111 class II). Previous studies 
of the young star population in the Trifid neb ula used 2MASS , 
ASCA, ROSAT ( Rho et al.l2001h and Chandra too et al.l2004l) . 

Finally, another possible contamination would be extragalac- 
tic sources. This effect is negligible based on the location of 
the Trifid nebula in the Galactic plane and close to the Galactic 
Center direction. 



3. Extinction law 

The general idea of the method to derive the extinction law is 
trivial. It requires the knowledge of the relative extinction in two 
wavelengths, A^lAy Aa2/Av, to evaluate the extinction at a third 
wavelength by measuring the slope of the star color distribution 
within a color-color plot m^i -m^2/'W/!2-'n^3- The color-color di- 
agram slope being equal to the color-excess ratio E;n-nl En-Ai 
the value of A^s/Av is inferred. Practically, this procedure is not 
necessarily that simple because of the scatter in the color-color 
diagram. The selection of a sub-sample as clean as possible is 
strongly recommended. Fig. [T^ is restricted to sources with an 
uncertainty smaller than 0. 15 mag at Ks, 3.6 and 4.5 yum although 
a large scatter is still present in the plot. In particular, the vertical 
elongation for [3.6] - [4.5] ~ and Ks - [3.6] > 2 is not real. 
It is due to a severe underestimation of the UKIDSS K^ uncer- 
tainty for some sources. Fortunately, the method presented in the 
next section makes their contribution negligible because of their 
random spatial distribution. 

Fig. [1] and several others in this paper represent plots as den- 
sity maps for clarity reason. Dark means high density. Unless 
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Fig. \.(a) Color-color diagram for stars with uncertainty smaller 
than 0.15 mag in the 3 bands, (b) Color-color diagram obtained 
using color maps instead of the source list. The exact same stars 
are used for both plots. The gray scale represents the density of 
the data points in the plot and the dots show the individual data 
points for low density. 



specified otherwise the overlaid symbols simply help to visual- 
ize the lowest density pixels that are hardly visible. 

3. 1. Color map making 
3.1.1. Method 

Although a slope could be measured in Fig. [T^, the scat- 
ter is large. The Trifid line-of-sight crosses various interstellar 
medium environments and all kind of stellar populations. This 
diversity, the high stellar density, and the relatively large dis- 
tance of M20 make our goals more difficult to reach than for 
nearby molecular clouds such as Chamaeleon, Taurus or Rho 
Ophiuchus. 

We propose a method based on the color maps that spectac- 
ularly improves the accuracy of the analysis. Fig. was built 
using the exact same sources as Fig. [T^ and both figures have 
the same axis, same scale. While each point in Fig. [1^ corre- 
sponds to a single source, points in Fig. [TJ) actually represent 
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Fig. 2. [3.6] - [4.5] color difference from sources de- 
tected simultaneously in two different sets of 3 filters, 
/:,(2MASS)/[3.6]/[4.5] and [3.6]/[4.5]/[5.8], versus color 



individual pixels taken from two color maps, one for each axis. 
Using the color allows us to keep the same axis as in Fig.[T^ but 
it is obviously equivalent to an extinction map assuming an ex- 
tinction law. The impressive difference between the two figures 
highlights the superiority of using maps rather than a source list 
directly. Using maps permits to average the color of several in- 
dividual nearby sources, and thus to reduce the scatter Different 
methods are available to build an extinction map. Star count 
methods a re sometirnes used on large scale an d for moderate 
extinction ('Cambresv"l999:'Dobashi et alJl2005 ^. This approach 
is definitely not adapted to a highly obscured distant cloud in the 
Galactic plane. Too sensitive to the foreground star contamina- 
tion the map would be noisy and seriously biased. Color-excess 
methods are better adapted. The first statistical method to map 
the exti nction from color excess, called NICE, was proposed 
by Lad a et al. d 19941). Improvem ent were proposed and called 
NICE R jLombardi & Alvesll2001h and then NICES T (Lombardi 
l2009h . A variant technique from ICambresv et al.l (i2002 ) pro- 
poses to replace the arithmetic mean color by the median color 
in cells and to fix the number of sources within a cell rather than 
its size. The latter step makes the method adaptive since the spa- 
tial resolution is always optimal, adjusted to the local density. A 
final convolution can be applied if a fixed spatial resolution is 
wished. This convolution minimizes a bias that would result in 
an underestimation of the extinction, compared to the usage of 
larger cells already at the final resolution size. Moreover the me- 
dian is a powerful tool to limit the effect of the residual contam- 
ination by foreground stars and young stellar objects as detailed 
in|Cambresv et al. (2002). 

All the color maps built for this work use the median color 
and follow the adaptive cell scheme with a fixed number of 3 
stars per cell. The variable cell spatial resolution being r,, a con- 
volution by a Gaussian kernel of -^1 



30 



^ 20 

X 

> 
< 

10 



1 


1 1 




' \ ' ' 'y-_ 
/ / - 


: 




case B 


/ X - 






/ 

■ / 


/ case A : 











































10 
Av(K3 



20 
[3.( 



30 



arcmin is applied to 
reach a final resolution of 1 arcmin over the whole map. 



Fig. 3. Extinction from H - color excess versus extinction 
from Ks - [3.6] color excess. By adjusting Ks - [3.6] to ob- 
tain a one-to-one regression line (solid line) we find A [3 f,] IAk, - 
0.6 11. The white iso-density contour delineates the region used 
to compute de regression line. The dashed line is the regression 
line for Av > 20 mag, it indicates the deviation at large extinc- 
tion. Case A and B refer to Table |2] 



tory to keep sources detected in all 3 bands to plot Fig. [T^, the 
Fig. [TJ) needs two color maps that can be built independently 
with only a 2 band detection criteria for each of them. Fig.|2]il- 
lustrates the color bias for [3.6] - [4.5] when it is obtained from 
two different sets of 3 band detection criteria, /r.5/[3.6]/[4.5] 
or [3.6]/[4.5]/[5.8]. To exaggerate the effect Ks is taken from 
2MASS which is not as deep as UKIDSS. The [3.6] - [4.5] aver- 
age color is affected by the detection criterion which makes 
the color bluer, and by the [5.8] detection criterion which tends 
to select redder objects. It is not a simple offset as shown in Fig.|2] 
because this effect is more important for redder color. In terms 
of extinction, taking sources detected in the 3 bluer wavelengths 
yields an underestimation of the extinction. Using deeper Ks data 
would obviously help since more sources detected at [3.6] and 
[4.5] would also be seen at Ks- However the effect cannot be 
avoided at large extinction because the Ks band suffer from twice 
the extinction than at [3.6]. 

It is therefore critical to compare only maps that are built 
with the exact same sample of stars. It implies the [3.6] - [4.5] 
color map needs to be built from a source list which contains 
stars detected in the 3 bands Ks, [3.6] and [4.5] if it is compared 
to a A'j - [3.6] map, or from a source list which contains stars 
detected at [3.6], [4.5] and [5.8] if it is compared to a [4.5] -[5.8] 
map. 

A total of 6 color maps are built for our study: H - Ks, 
Ks - [3.6], [3.6] - [4.5], [4.5] - [5.8]. Each map is obtained with 
sources detected in 3 bands at the same time. There are 2 ver- 
sions for the Ks - [3.6] map and the [3.6] - [4.5] map because 
the third band can either be the bluer or the redder one. 



3.1.2. Selection bias 

Once a clean star sample is produced, i.e. filtered from fore- 
ground stars, YSOs and noisy photometry sources, the color 
maps are still subject to a selection bias. While it is manda- 



3.2. Results 

Although the extinction law has been for long considered as uni- 
versal in the infrared it is now admitted there are variations at 
wavelength longer than 2 yum. The spectral range that is pseudo- 
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universal is actually restricted to 1-2 /im where the extinction 
is independent or very little dependent of the environment, i.e. 
dust temperature an d density. Consequently, An/Ay - 0.175 
and AkJAv = 0.1 12 dRieke & Lebofskvl 19851) are assumed con- 
stant. It corresponds to j3 - 1 .7 in the expression A^ oc A^^ that 
is often used to describe the near-infrared extinction law. Values 
from 1.6 to 1.8 being reasonable, the effect of such a systematic 
uncertainty is discussed further in this section. The H - map 
can be converted into an extinction map using 



(H - K,) - {H - K,)o 
Ah/Av -AkJAv 



= 15.87 Eh-k, 



The extinction law is normalized to Ay mainly for historical 
reason but it is more convenient to normalize to Kg. The re- 
lation between the extinction and the color excess becomes 
Ak^ - 1.78 Efj-K,- The reference color (H - Ks)o - 0.54 mag is 
estimated from the area with the smallest values within the map, 
about 50 arcmin south-west from the densest core of the cloud. 
This is much redder than the intrinsic stellar color because it 
includes the diffuse extinction along the line-of-sight. The ref- 
erence color is relevant to provide an absolute extinction for the 
Trifid alone and it is discussed in the Sect. 14.21 It has absolutely 
no effect on the extinction law estimation nor on the relative ex- 
tinction comparison from different parts of the field. 

Fig.[3]represents Ay from H - versus Ay from - [3.6]. 
Since Ay - AkJQ.112 is assumed, and by definition Ak, - 
Ek,-[3.6]/{^ - A[3.6]/Ajf ), the only free parameter to adjust to ob- 
tain a slope of 1 is the value of Ap ej/A/f . After the adjustment, 
it yields the value of Ap ej/A^-, = 0.611. Fig.|3]is presented as a 
density map for clarity and also because the density information 
is taken into account in the fitting process to exclude the low den- 
sity pixels. Basically the minimum density level is fixed to 20% 
of the maximum density in the map. The color-coding used in 
the density plots is in log-scale so the 20% always corresponds 
to the dark region. The white iso-density contour delineates the 
region selected to compute the regression line in Fig|3] Any out- 
lier is therefore fully eliminated with this technique. 

Once A[3.6]/Ajf^ is determined, A[4.g]/A;i:,, is evaluated with a 
plot of Ay from - [3.6] versus Ay from [3.6] - [4.5]. A third 
similar plot will then provide Ap gj/A/f^ (see Fig.|4]i. The results 
areA[4.5]/A;f, = 0.500 and Ap.gj/A/r,, = 0.398. 

Since the pixels falling in the low density regions of the Ay 
vs. Ay plot are ignored the fitting is performed on regions where 
the visual extinction is roughly smaller than 15 mag. At large 
extinction a striking deviation of the slopes is observed and out- 
lined by the dashed lines in Figs.[3]and|4] The changes of slope 
occur at Ay 20 mag (Ak, ~ 2.2) and suggests a transition 
phase in the dust grain population. At such a high column density 
this result is unexpected. Contrary to studies in the far-infrared 
wavelengths which rely on the dust emission, there is no temper- 
ature degeneracy issue since it is based on stellar color mapping. 
A drawback of this technique is its limitation to the near-infrared 
spectral range. At wavelengths longer than ~5 //m the number of 
detected stars is too low to have an acceptable spatial sampling 
of the molecular clouds while shorter wavelengths are so sensi- 
tive to the extinction that the cloud dense regions become totally 
opaque. A direct consequence on this limitation is that no ev- 
idence for a slope change can be seen in a J - H vs. H - Ks 
diagram because the visual extinction cannot go deeper than 
~20 mag with J-H. That prevents any conclusion on the extinc- 
tion law variation at shorter wavelengths than 3.6 fim. Table |2] 
summarizes the results obtained from the slope measured in both 
Figs. [3] and |4] We consider two different cases for Ap e] /A /j-^ at 
Ay > 20 mag: 
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Fig. 5. Extinction law. The two red lines corresponds to the cases 
described in Table|2]for Ay > 20 mag. The shaded area indicates 
the systematic error, for case B only, resulting from Ah /Ak, un- 
certainty which can range from 1.52 (J3 - 1.6) to 1.60 (J3 = 1.8). 
At longer wavelengths than Ks, /3 = 1. 6 gives larger A^/A^.. 
Diamonds are for Ay < 15 mag. The IWeingartner & Draind 
(2001) model is plotted in gray with dashed and plain lines for 
Ry - 3.1 and Ry = 5.5, respectively. The axis scale is linear for 
X and logarithmic for y. 



- case A: A[3 6]/Aa:, - 0.611 is independent from the extinc- 
tion, the value remains identical than at lower extinction 
(solid line in Fig. [3). 

- case B: A[3 6]/Aa:^ - 0.678 which is the value actually ob- 
tained from the slope of the dashed line in Fig. [3] 

The statistical errors on the slope estimations are small and 
so are the errors on the computed A^/Ak^- The systematic er- 
ror on Ah/Ak, is the major source of uncertainty. We assume 
a value of 1.563 which corresponds to a power-law index of 
P - l.l for the extinction law. Our results from Table |2] are 
plotted in Fig.|5] The shaded area delineates the systematic un- 
certa inty for P varying fro m 1.6 to 1.8. Using (3 - 1.6, such 
as in Chapman et al. (2009), increases AjJAjr to the upper bor- 
der of the shaded area. Nishiy ama et a 11(120061) reported a higher 
valu e of ;8 = 1 .99 which is not com patible with our observations. 
The IWeingartner & Drain3 (12001 ) (hereafter WDOl) model is 
overlaid for comparison. The variations in the extinction law 
found in this work are rather small compared to the dust evo- 
lution from Ry - 3.1 to Ry - 5.5. The variation being observed 
at visual extinction larger than 20 mag is likely the main rea- 
son for this dust transition to be unknown. In case A the ex- 
tinction parameters at 4.5 and 5.8 jum are smaller at large ex- 
tinction than at small extinction which is not realistic accord- 
ing the WDOl model. This case teaches us the small deviation 
observed in Fig. [3] cannot be ignored as it would make longer 
wavelength results inconsistent. Case B is more likely even if 
it yields a larger value than reported in the literature until now. 
The largest value prior to this work was from IChapman et al.l 
(2009) with A[3.6]/Aa:, = 0.64. H owever their numbers are over- 
estimated compared to ours or to lFlahertv et alj (l2007l) because 
they chose a smaller yS. Their results come from an average for 
Ophiuchus, Perseus and the Serpens molecular clouds which 
all form low-mass stars whereas the cloud assoc iated with the 
Trifid is clearly more massive and forms OB stars. iFlahertv et aP 
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Fig. 4. Each plot represents the extinction correlation from two different color excesses. The two parameters A [4 5] /A = 0.500 (left 
plot) and A[5 8]/A/i-^ = 0.398 (right plot) are adjusted to obtain a one-to-one regression line (solid line). The dashed lines trace the 
deviation at large extinction, the corresponding slopes in these plots are 0.55 + 0.02 and 2.52 + 0.09, respectively. 

Table 2. Color excess ratios and resulting extinction law parameters. The two rows for Ay > 20 mag represent different cases 
for A [3 6] /A jf^: case A propagates the same value as for small extinction, case B takes the actual measured value. We assume the 
power-law index /? = 1.7 (i.e. Ah/Ak, - 1.563) and results for y8 = 1.6 and 1.8 are provided for case B. The main results from the 
literature for dense interstellar medium directions and from a dust model are indicated for comparison. 





Eh-K, /Ek, -[3.6] 


Ek, -13.6]/ E[3_(,]-[4.5] 


-£[3.6]-[4.5]/£'[4.51-[5.81 


Ah/Ak, 


A]3.6]/Ak, 


A]4.5]/Ak, 


A[5.8]/Ajf, 


Av < 15 


1.446 + 0.002 


3.508 ±0.011 


1.087 ±0.007 


1.563 


0.611 ±0.002 


0.500 ± 0.003 


0.398 ± 0.005 


Av > 20 
















case A 


1.446 ±0.002 


1.915 ±0.018 


1.496 ±0.089 


1.563 


0.611 ±0.002 


0.408 ± 0.003 


0.272 ±0.011 


case B 


1.749 ±0.003 


2.317 ±0.027 


1.810 ±0.091 


1.563 


0.678 ± 0.006 


0.540 ± 0.003 


0.463 ± 0.007 


case B, yS = 1.6 








1.520 


0.703 ± 0.005 


0.573 ± 0.003 


0.500 ± 0.007 


case B,/3 = 1.8 








1.601 


0.656 ± 0.006 


0.507 ± 0.003 


0.424 ± 0.007 



Roman-Zuniga et al. (2007) for Barnard 59 1.55 0.619 ± 0.029 0.512 ±0.022 0.459 ± 0.024 

Flaherty at al. (2007) for Orion 1.55 0.636 ± 0.003 0.540 ± 0.003 0.504 ± 0.003 

Chapman et al. (2009) for Ophiuchus, Perseus and Serpens 1.52 0.640 ± 0.030 0.530 ± 0.030 0.460 ± 0.030 

Weingartner &Draine (2001) with = 3.1 0.40 0.25 0.17 
Weingartner &Draine (2001) with Ry = 5.5 060 049 0.40 



(l2007h investigate the high-mass star-forming region of Orion 
and obtain a similar number of 0.636 + 0.003, significantly be- 
low our A[3 6]/Aa:, - 0.678 + 0.006. They did find smaller values 
at all wavelengths in low-mass star-forming regions. The pres- 
ence of massive star, in addition to a high density, may therefore 
be relevant. 

WDOl attribute the variation of the extinction law with Ry 
to grain growth only, although ice mantles onto dust grains con- 
tribute to the optical property variations. The main ice features 
relevant to our study originate from H2O and CO2 ices with ab- 
sorption at 3.0 jum and 4.27 fj.m, respectively. Both H2O and CO2 
ices are detected at visual extinction below 4 mag ( Whittet et alJ 
[2001 ; Bergin et al. 2005) where the WDOl model is still valid. 
This indicates ice is not sufficient to significantly modify the ex- 
tinction law in the Spitzer/IRAC filters. We can also rule out 
the possibility that ice features saturate at Ay «j 20 mag be- 
cause they are observed at larger extinction without any satura- 
tion (A. Boogert, pri vate communication) . An example of spec- 
trum is presented in IWhittet et al.l ( Il998l Fig.l) for a star with 
Ay - 21 mag. It worth noting the H2O ice band at 3.0 yum 



peaks outside the [3.6] IRAC filter (3.2-3.9 fim). Beside the ab- 
sorption features, ice mantles make the grains more sti cky, help- 
ing the coagulation process. Using Spitzer/IRS spectra lMcClurd 
(2009) found evidence for a flattening of the extinction law at 
Afc > 1 mag which would be due to the coagulation of ice- 
mantled grains. This suggests ice makes the coagulation process 
more efficient but it is not the dominant parameter since the co- 
agulation starts well after ice-mantles appear. The abrupt change 
in the extinction law at Ay x 20 is not understood and we can 
only speculate that the high density in the Trifid cloud might 
allow a new structuring of the dust grains or that the shocks ob- 
served in the Trifid cloud could be determinant. In a recent work 
Boogert et al. (submitted to ApJ) also reported a flatter extinc- 
tion law than WDOl based on a spectrum from 3 to 24 fj.m for a 
star with 25 mag of visual extinction. 
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Fig. 6. Left: large scale extinction map of the Trifid region. Right: GLIMPSE 8 jim image of the Trifid nebula with iso-extinction 
contours. In both figures contours correspond to Ay - [10, 15, 20, 40, 60] mag, the dashed line follows the Galactic plane. The red 
shaded regions delineate the 5 1.25 > 10 mJy/H" beam Trifid Condensations (TC. ICernicharo et alj|1998h studied in Sect. 14.31 at a 
1 arcmin resolution. 



4. From color to extinction 

4.1. Composite map 

The extinction law parameters derived in the previous section 
allow us to convert the color maps into extinction. To proceed we 
follow the more realistic case B for A^/A^r^ when Ay > 20 mag 
(see Tabled and a reference color of H - Kg = 0.54 mag. The 
associated zero points derived for the other colors are Ks-[3.6] - 
0.497, [3.6] - [4.5] = -0.017 and [4.5] - [5.8] = 0.187 mag. 
The sensitivity of each map relies on the wavelengths used to 
build the color maps. Shorter wavelengths are more sensitive too 
smaller extinctions. The final extinction map presented in Fig.|6] 
is a combination defined as follow: 

( K, - [3.6] if Av < 15 mag 
Ay from \ [3.6] - [4.5] if 10 < Ay < 80 mag 
[ [4.5] - [5.8] if Av > 60 mag 

The transition between two maps is progressive and linear within 
the overlapping extinction ranges. A maximum visual extinc- 
tion of ~80 mag is reached. In complement to the extinction 
map Fig. |6] shows a close-up on the Trifid nebula observed at 
8 fim with the iso-extinction contours superimposed. The parent 
molecular cloud is an elongated structure with its densest region 
lying just West of the famous optical nebula. 

4.2. Tridimensional analysis, distance estimation 

The extinction estimation relies on background stars. The mean- 
ing of background is not that trivial in the direction of the Trifid 
nebula since the line-of-sight crosses the whole galactic disk at 
only 1 kpc of the Galactic center The complexity of the line-of- 
sight and the distance of the stars used for the analysis are criti- 
cal and can change our interpretation of the color maps. Even the 
distance of the nebula is still uncertain. To address these issues 



we performed a tridimensional analysis which aims to separate 
the different components along the line-of-sight. 

A three dimensional stellar model of the Galaxy ("Robi n et"al] 

1 20031) is used to simulate the stellar distribution along the line- 
of-sight towards the Trifid nebula. The simulated stars are not 
subjected to extinction and so they provide an estimate of 
the intrinsic color and distance of the stars along this line-of- 
sight. By comparing the color distribution of the unreddened, 
simulated, stars to the observed stellar color distribution it is 
possi ble to infer the extin ction distribution along the line-of- 
sight ([Marshall et alJl2006l) . We use observations from 2MASS, 
UKIDSS and GLIMPSE and use the Ks-[3.6] color to probe the 
line-of-sight extinction. This technique is relatively independent 
of the stellar counts and reflects the changing stellar color as a 
function of distance (see [Marshall et al1l2006l for more details). 
Furthermore there is no need to remove foreground stars as we 
are comparing the observations to a model where all stars are in- 
cluded. However, care needs to be taken to ensure that the model 
is indeed representative of the observations, so the simulated cat- 
alog is cut at the completeness magnitud e of the observa t ions in 
the Ks, [3.6] and [4.5] bands. Recently, ' MarshaU et all (|2009|) 
improved on this method using a genetic algorithm to deduce 
the line-of-sight extinction necessary to reproduce the observed 
stellar color distribution along a given line-of-sight. This version 
is the one we chose to use as it provides higher angular resolu- 
tion (~3 arcmin compared to 15 arcmin previously) and more 
sensitivity to extinction features within 2 kpc from the Sun. The 
resulting distances are subject to an uncertainty of roughly 15% 
but this may be higher at large distances and through high col- 
umn densities. This results in a smearing out of extinction fea- 
tures along the line of sight, and could explain why, in Fig.|2l a 
narrow extinction peak is found at 6.5 kpc when the extinction is 
low but this same peak is found to be wider through the higher 
density lines-of-sight crossing the Trifid nebula. At such a dis- 
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Fig. 7. Extinction (top) and star density (bottom) along the Trifid 
line-of-sight. The black line is the average value for the region 
where Ay > 20 mag in Fig. |6j the gray line is for the region 
where Ay e [0, 1] mag. The dashed lines indicate the Trifid lo- 
cation at a distance of 2.7 + 0.5 kpc. 

tance the interstellar material can be related to several Galactic 
structures which are the molecular ring (about 3-4 kpc from the 
Galactic center), the start of the Norma Arm, and the extremity 
of the Galactic bar 

The results of the 3D modeling is presented in Fig. |7] The 
extinction profile for the region where Ay > 20 mag (see Fig.|6]l 
provides us with a direct estimation of the cloud distance. Fitting 
a Gaussian yields a distance of 2.7 + 0.5 kpc. The extension of 
the cloud along the line-of-sight is likely from 50 to 100 pc so 
the Gaussian standard deviation represents the uncertainty inher- 
ent to the 3D m odeling, not the c loud size. The literature often 
refers to 1 .7 kpc (^ Lvnds et aljl98 5) for the Trifid nebula distance 
althou gh a serious alternative is proposed by Kohoutek et al] 
(Il999h with a distance between 2.5 and 2.8 kpc. Both estima- 
tions are based on the photometry of the central stars. The sec- 
ond work is more complete as it uses more stars and includes 
spectroscopy and accurate photometry. Our analysis of the tridi- 
mensional stellar distribution unambiguously supports the long 
distance hypothesis which place the Trifid in the Scutum Arm 
(Iwii3l2008.) rather than the Sagittarius Arm as generally re- 
ported. In the following we assume d = 2.7 kpc. 

For such a distance, the foreground star number is expected 
to be at least 37 x 10^ deg"^ (see Fig. [ill, or 11 x 10^^ deg"^ if 
we restrict the sample to 2MASS sources. This is in good agree- 
ment with the analysis presented in Sect. 12.21 where we found a 



total foreground star density of at least 37 x 10^ deg"^ towards 
the dense region of the cloud and in agreement with the distribu- 
tion of 8600 stars deg"^ detected by 2MASS which was found to 
be uncorrected with the extinction, i.e. foreground to the cloud. 
It worth noting the Sect. 12.21 results are also incompatible with 
the short distance hypothesis. The '^CO velocity obtained from 
iDame et 30 (12001) shows a velocity structure consistent with our 
extinction distribution which corroborates our estimation. The 
Trifid molecular clouds LSR ve locity is about 18-20 km s "^. 
Assuming a Sun kinematic from lMcMillan & BinnevI (1201 Ol) . a 
distance to the Galactic center of 7.8 kpc and a rotation speed 
of 247 km s"' the Trifid velocity corresponds to a kinematic dis- 
tance of 3.2-3.4 kpc which is even larger than our estimation 
based on the stellar population model. 

The 3D analysis estimates the diffuse extinction in this direc- 
tion to ~3 mag kpc"' and it shows that no other dense cloud sig- 
nificantly contaminates the Trifid line-of-sight. The stellar frac- 
tion at distance larger 8 kpc represents only 7.7% of the back- 
ground stars at large extinction and 2.6% at small extinction. In 
both cases their median distance is ~5.8 kpc. The zero point cho- 
sen for the extinction map (Fig.|6l) is validated by the 3D study 
as 1) there is no peak at ~3 kpc at low extinction and 2) it com- 
pensates for the diffuse extinction on the line-of-sight and the 
contamination at 6.5 kpc. 

4.3. Cloud mass 

The cloud mass is derived from the extinction map using the 
following relation (Dickman 1978.) : 

M^iacD^fi^y Av(i) (1) 

Ay ^-j-' 

where a is the angular size of a pixel map, d the distance to 
the cloud, fi the mean molecular weight corrected for the he- 
lium abundance, and / re presents a pixel map. W ith the dust- 
to-gas ratio proposed bv ISavage & MathisI (Il979l) . Nh/Av - 
1.87 X 10^' cm"^ mag"' (Nh - Nhi + Nh^) and a distance of 
2.7 kpc, we obtain a mass of 5.8 X 10^ Mq. This value corre- 
sponds to twice the Orion (A+B) clouds mass (ICambres'vlll999l: 
IWilson et alJ2005h . Figure [8] shows how the integrated mass dis- 
tribution varies with the extinction within the cloud. Regions 
with visual extinction greater than 20 mag account for less than 
4% of the total mass of the cloud. The index of the cumulative 
mass distribution, a, changes at Ay ~ 25 mag by a factor of 2. 
For M(Ay) being the mass enclosed in the contour of extinction 
Ay, we have: 

MiAy) = M(0) X lO"""-*" (2) 

. f Q- = 8.39 X 10-2 forA^e [8,20] 
^™ I a = 3.76 X 10-2 for Ay e [30,60] 

The variation of a is not a side-effect of the map combina- 
tion technique nor the application of a corrected extinction law 
at large Ay since the variation is also observed for a raw ex- 
tinction map obtained only from the [3.6]-[4.5] color-excess and 
assuming a constant extinction law. It is however difficult to in- 
terpret the index variation as a consequence of the dust transition 
responsible of the extinction law evolution. The two events oc- 
curring at the same level may just be a coincidence. The Nn/Ay 
value is generally assumed to be universal although some au- 
thors su ggested a d ecrease by a factor of 3 in deeply embedded 
regions ( IWinston"et al. 2010). Dividing the mass estimate by any 
factor above a given extinction threshold would only produce an 
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Fig. 8. Cumulative mass distribution for d-2J kpc. The a index 
is defined in Eq.|2] 



offset in Fig.|8] not a change of slope. A more complex evolution 
of Nh/Av with the extinction would affect differently the Fig. [8] 
distribution but it seems unlikely to only produce a sharp slope 
change. 

The millimeter wavelength emissio n from th e cold dust 
continuum has been mapped by Cernicl wo et al.l (Il998h and 
iLefloch et al.l (l2008l) . They have detected several dust conden- 
sations from the 1 .25 mm emission. For simplicity we designate 
these dust condensations TC NN as in the original papers. This is 
fine is the Trifid context but obviously ambiguous in the general 
context and the full SIMBAD name is actually [CLC98] TC NN. 
For the more extended condensations we compare the mass from 
the dust infrared extinction and 1.25 mm emission. The results 
are presented in Table|3] The mass from extinction is obtained by 
integrating the area defined in the map at 1.25 mm, convolved at 
1 arcmin resolution, by the contour at 5 1.25 - 10 mJy/H" beam 
(red shaded area in Fig. |6]l. Because of the observational tech- 
nique in the millimeter wavelengths, a spatial filtering removes 
the structure at scale larger than 4.5 arcmin in the dust contin- 
uum emission. To compare this map with the extinction map we 
scale it by adding the median extinction measured in the vicinity 
of each TC region. Basically the median extinction is measured 
in the area between the contours at 5 1.25 - 5 and 10 mJy/H" 
around each condensation. The resulting offsets are given in the 
column labeled outer of Table [3] The mass from the map at 
1.25 mm is derived assuming a dust temperature of 15 K and a 
dust opacity of a:i 25 = 4.08 x 10"' cm^ gr"' following the model 
B for Rv - 5.5 by WDOl. The corresponding extinction cross 
section per hydrogen atom is C1.25 = 7.336x10"^^ cm^ H"' after 
correction fo r the helium abund a nce and assum ing a gas-to-dust 
ratio of 125 jLi & Drainj|2001h . lLefloch et al.l (|200 8) measured 
temperatures from 20 to 30 K using CO observations. The gas ki- 
netic temperature and the dust temperature are however different 
in such regions and an arbitrary Tdust - 15 K is more appropriate. 
The masses we derived from the dust emission are systematically 
smaller. Two parameters, the temperature and the emissivity, can 
be adjusted to reach an agreement with the masses from extinc- 
tion. While TCOO and TO only need a small increase of the dust 
emissivity at 1.25 mm, or to decrease the temperature from 15 to 
13 K, TC3/4/5 require a more significant emissivity increase by 
a factor 1.7 to 2.7. Using only the temperature as free parameter 



would lead to 11 to 8 K. Such low temperatures have already 
been reported by Desert et"an (l2008h from submillimeter point 
sources observed by the Archeops experiment. However, the 
condensations we are dealing with cover an extended field from 
3 to 12 arcmin^. Assuming a minimum temperature of 12 K for 
TC3/4, an emissivity enhancement by a factor 1 .9 would make 
the mass estimates in agreement. TCOO and TCO are starless con- 
densations while TC3/4/5 form stars and undergo shocks. Our 
analysis suggest the dust emissivity at millimeter wavelength 
is enhanced for these condensations, especially for the densest 
ones. Complementary observations at longer wavelengths, from 
Herschel and Planck, will permit to better understand the physi- 
cal properties of the interstellar medium in the Trifid region. We 
plan to follow up this study as soon as the data from the Hersche l 
Infrared Galactic Plane Survey (Hi-G AL. iMolinari etal.ll2010bh 
are available. 

5. Conclusion 

We performed a detailed analysis of the extinction in the molecu- 
lar cloud associated with the Trifid nebula. The direction and the 
distance of this cloud is a major difficulty compared to the other 
similar investigations published in the literature. To overcome 
this specificity we proposed an original method which combines 
the color-excess mapping with the color-color diagram analysis. 
The gain in sensitivity with this method allowed us to measure 
to extinction law at 3.6, 4.5 and 5.8 fim and to unambiguously 
detect a variation at large absorption through an abrupt change 
of slope in the color-color plots. We interpreted this result as an 
evidence for a rapid dust evolution at high density like a new 
dust transition phase in the dense interstellar medium. 

The extinction law parameters, A^/A^^, are found to be larger 
in the very dense cores of the cloud, in agreement with dust 
models. Our values for Ay > 20 mag are however larger than 
predicted by WDOl for dense regions (Ry = 5.5) and larger than 
previously reported in the literature at 3.6 yum. For Ay < 15 mag, 
which does still refer to dense material, our values remark- 
ably match those predicted by the model for Ry - 5.5, as ex- 
pected. Our study is not sensitive to the diffuse medium for 
which Ry = 3.1 with a typ i cal vis ual extinction as low as 1 mag. 
Moreover, IZasowski et al.l (l2009l) found a correlation of the ex- 
tinction law for the diffuse interstellar medium with the galac- 
tocentric radius. They concluded the extinction curve follows 
the dust model with Ry - 5.5 for a galactocentric radius con- 
sistent with the Trifid distance. It suggests the so-called diffuse 
medium already contains big grains in this part of the Galaxy. If 
the diffuse medium at the Trifid distance behaves like the dense 
medium in the solar neighborhood we can imagine the Trifid 
dense cores are also affected by the distance with a flatter ex- 
tinction law. 

Using our varying extinction law we built a composite ex- 
tinction map of the Trifid nebula region. The maximum visual 
extinction is 80 mag at 1 arcmin resolution. We have also inves- 
tigated the matter distribution along the line-of-sight since sev- 
eral small clouds could have mimicked the presence of a single 
giant massive dark cloud. We found no significant contribution 
of other clouds in this direction and we were able to estimate the 
distance of the Trifid nebula at 2.7 + 0.5 kpc. 

The Trifid molecular cloud is about twice the mass of the 
Orion molecular cloud. The comparison of the dust extinction 
and 1 .25 mm continuum emission in the densest cores suggests 
the dust emissivity is probably enhanced by a factor of 2-3 in 
these regions. This cloud appears as a young and massive ver- 
sion of the Orion or Rosette regions. The Trifid represents an 
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Table 3. Properties of the more extended Trifid cores detected at millimeter wavelength (iLefloch et alJ l200"8h . Masses are 
expressed assuming a distance of 2.7 kpc. The mass from dust continuum emission is derived using Tdust = 15 K and 
/ci.25 - 4.08 X 10"' cm^ gr"' (i.e. C1.25 - 7.336 x 10"^^ cm^ H"'). By modifying these parameters the mass can be adjusted 
to the mass estimated from the dust extinction. The two last columns give these parameter values, providing only one of the two 
parameters is modified at the same time. 



Name 
[CLC98] 


Mass^ 
(M~ 


Mass^ 
(M~ 


peak 


Av (mag) 
median 


outer 


Surface 
(arcmin^) 


K 
(X) 


(K) 


TCOO 


936 


911 


27 


22.7 


20.0 


3.5 


1.3 


13 


TCO 


1607 


1570 


25 


20.3 


16.9 


6.9 


1.2 


13 


TC3 


3535 


2783 


78 


40.4 


27.6 


7.3 


2.6 


8 


TC4 


1149 


945 


52 


32.8 


23.7 


3.0 


2.7 


8 


TC5 


3543 


3087 


53 


25.5 


18.1 


11.6 


1.7 


11 



" Mass from the dust near-infrared extinction. 

' Mass from the dust continuum millimeter emission. 



earlier evolutionary status in the star-formation process with a 
first generation of OB stars but no significant embedded clus- 
ter yet. Several protostars in TC3/4 are likely the precursors of 
such clusters. Longer wavelengths observations from Herschel 
and Planck are needed to better understand the dust properties 
and the interaction between the interstellar medium and the OB 
stars. 
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